esp@cenet document view 



Page 1 of 1 



View-to-view image correction for object motion 



Patent number: 
Publication date: 
Inventor: 
Applicant: 
Classification: 

- international: 

- european: 
Application number: 
Priority number(s): 



FR2602602 
1988-02-12 

NOWAK DAVID JAMES 
GEN ELECTRIC (US) 

G06F15/66; G06F15/42; A61B6/03; G01T1/166 
G01T1/29D4; G06T11/00T 
FR1 987001 0784 19870730 
US1 986089531 6 19860811 



Also published as: 

US4858128 (A 
JP63088688 (A 
GB21 94870 (A 
DE3726595 (A 



Report a data error he 



Abstract not available for FR2602602 

Abstract of corresponding document: US4858128 

An imaging system of a type producing a plurality 
of X-Y matrices representing projection views of 
an object for later reconstruction using parts of a 
plurality of the views employs correlation 
techniques for determining the amount of motion 
along X and Y axes. The detected motion may be 
removed by suitably shifting later data to align it 
with earlier data, or vice versa. In one 
embodiment, all brightness values parallel to the 
X axes of two views are summed to develop two 
one-dimensional arrays containing information 
about the brightness patterns in the Y direction. 
One of the two one-dimensional arrays is shifted, 
and the two one-dimensional arrays multiplied, 
pixel-by-pixel, and the products are summed to 
provide a measure of similarity between the 
brightness patterns in the one-dimensional 
arrays. Repeated shifts produce a correlation 
function the shift position of whose maximum is a 
good approximation of the amount of motion 
occurring between the two views. The detected 
motion may be used to alert the operator that the 
data is faulty or to permit automatic alignment of 
all data from a scan. An identical process may be 
performed along the Y axis to detect and correct 
motion along the X axis, in certain restricted 
cases. 
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Dispositif et procede de correction d'image vue a vue pour le mouvement d'un objet. 
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Un systeme d'imagerie produisant une multitude d'ensem- 
Bles X-Y representant des vues en projection d'un objet pour 
reconstruction ulterieure en utilisant des parties des vues 
emploie des techniques de correlation pour determiner le 
mouvement suivant les axes X et Y. Le mouvement detects 
peut etre elimine en proc&Jant au decalage 58 d'une donnee 
ulterieure pour f aligner a une donnee anterieure, ou vice versa. 
Dans un mode de realisation, toutes les vateurs de brillance 
paralleles a 1'axe X des deux vues sont addrtionnees pour 
developper deux ensembles a une dimension contenant une 
information sur les diagrammes de brillance dans la direction 
Y. L'un des deux diagrammes est decale, et les deux ensem- 
bles sont multiplies, element par element et les produits 
addttionnes pour foumir une mesure de la similarity entre les 
diagrammes. Des decatages repetes produisent une fonction 
de correlation, la position de decalage du maximum etant une 
bonne approximation du mouvement produit entre les deux 
vues. On peut utiliser les mouvements. detectes pour signaler a 
I'operateur que la donnee est defectueuse ou pour permettre 



un alignement automatique de toutes les donnees d'un ba- 
layage. On peut executer un processus identique pour I'axe Y. 

Application aux systemes dlmagerie, par exemple a llmage- 
rie par tomographie informatique aux rayons-X. 
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La presente invention concerne I'imagerie et, plus 
particulierement, la detection du mouvement d'un objet, et 
sa correction, dans un systeme d* imager ie ou des vues suc- 
cessives sont prises a partir de positions differentes, 

5 chacune representant une projection a deux dimensions de 
1' objet a trois diraensions- 

Bien qu'on puisse employer la presente invention 
dans d'autres types de systeraes d* imager ie tels que, par 
exemple dans 1 1 imager ie radiographique a tomographie inf or- 

10 matique, le memoire suivant est relatif au cas d'un systeme 
tomographique a emission dans le but d'en rendre la descrip- 
tion plus concrete . 

Un type de systeme tomographique a emission pre- 
sentant de l'interet est le systeme de tomographie informa- 

15 tique a une seule emission de photons (TISEP), dans lequel 
un emetteur de rayons gamma de faible niveau est injecte 
dans le corps d'un patient. L' emetteur classique de rayons- 
gamma est d'un type qui se propage pref erentiellement vers 
un organe dont on doit produire 1* image. Un detecteur de 

20 rayons-gamma, plan, de grande surface, detecte les rayons- 
gamma emis par le corps perpenticulairement a son plan. 
Cette information est stockee numeriquement sous forme d'une 
image dans un reseau M x N d* elements appeles pixels ou 
elements d* image. Les valeurs de M et N sont convent ionnel- 
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lement egales l'une a X' autre, et sont couramment de 64 ou 
128 unites, ou pixels, dans les deux dimensions de 1' image. 

Un sy steme TISEP emploie une multitude de vues, 
chacune etant prise en placant un detecteur parallelement a 

5 un axe fixe, et suivant un certain angle par rapport a cet 
axe. Zb 1 angle est incremente par etapes egales entre les 
vues. Les images ainsi captees sont traitees sur ordinateur 
afin de construire des images de tranches transaxiales du 
corps. Le traitement informatique utilise des parties de 

10 chacune des vues se succedant pour reconstruire chaque 
tranche transaxiale. Si le corps dont on fait une image 
change de position entre des vues successives, la donnee due 
a une vue n'est pas en alignement avec les donnees provenant 
d'autres vues. II en resulte que les images reconstruites d 

15 partir des donnees peuvent etre floues, ou contenir des 
elements de distorsion dans le processus d'imagerie qui ne 
representent pas les caracter istiques reelles de 1'objet 
examine. 

Afin de reduire au minimum le dosage du rayonne- 
20 ment auquel le patient est expose, les materiaux injectes 
qui emettent des rayons-gamma ont une radioactivite relati- 
vement faible. II en resulte que la production de chaque vue 
necessite un temps important, tel que, par exemple, 40 
secondes environ. Si l'on souhaite prendre un total de 64 
25 vues sur un arc de 360°, celles-ci etant espacees d* environ 
5,6°, il faut alors environ 40 minutes pour completer le 
processus complet d'imagerie. Un flou ou une distorsion peut 
se produire lorsque le corps examine se deplace d'une dis- 
tance de I'ordre d'un element d 1 image. Un element d' image 

2 

30 typique a une surface d'environ 1/2 • cm . L' immobilisation 
du corps humain sur une surface de moins d'un demi centi- 
metre carre pendant une duree de 40 minutes est difficile a 
obtenir, voire impossible. C'est pourquoi il est courant 
qu'on assiste a un mouvement du corps et a la degration qui 

35 en resulte pour les images. 



On peut detecter l 1 existence du mouvement du corps 
a partir de Xa donnee enregistree dans une procedure dite a 
enregistrement en mode vertical, dans laquelle on affiche 
1" ensemble des vues dues au balayage, l'une apres 1' autre, 
afin de creer une image simulee de cinema a partir de la- 
quelle l'operateur peut observer s'il s'est produit une 
valeur inacceptable du mouvement du corps pendant la col- 
lecte des donnees. Cette technique est f ondamentalement une 
technique de controle de qualite permettant de determiner si 
la donnee recueillie est utilisable. Si la donnee n'est pas 
u tilisable, 1 'alternative est soit de decompter de fa<jon 
appropriee les donnees d' image recueillies pendant 1 1 inter- 
pretation des images soit de repeter la collecte des don- 
nees. Cette technique n'apporte pas un moyen permettant de 
corriger les donnees pour eliminer les erreurs dues au 
m ouvement. De plus, la determination qu'effectue l'operateur 
est au moins en partie subjective et done encline aux er- 
reurs. 

Une autre facon de detecter le mouvement du corps 
est appelee sinogramme, a savoir une image creee en affi- 
chant la donnee recueillie qu'on utilise ulterieurement pour 
construire une tranche transaxiale. L'operateur se base sur 
1 'observation des elements faconnes du mouvement du corps en 
detectant visuellement certaines distorsions dans 1* image du 
sinogramme. Comme avec 1 1 enregistrement en mode vertical, il 
s'agit la avant tout d'une technique de controle de qualite, 
qui ne permet pas la correction des erreurs dues au mouve- 
ment. Celle-ci souffre similairement de 1 'obligation de 
proceder a des jugements subjectifs de la part de l'opera- 
teur. 

Une autre technique, telle que celle decrite, par 
exemple, dans un article de J.S. Fleming intitulee "A Tech- 
nique for Motion Correction in Dynamic Scintigraphy" (Tech- 
nique de correction du mouvement dans la Scintigraphic 
dynamique) paru dans la revue European Journal of Nuclear 
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Medicine (1984) volume 9, pages 397-402, fait appel a des 
sources ponctuelles d* emission de rayons-gamma qui sont 
appliques au corps dont on fait 1* imager ie. Les sources 
ponctuelles font l'objet d'une imagerie en meme temps que le 

5 reste du corps. I*e mouvement detecte des sources ponctuelles 
peut etre utilise dans une" procedure de controle de qualite 
et fournir des donnees suffisantes pour appliquer des fac- 
teurs de correction manuelle a certaines des donnees affec- 
tees- Cette technique souffre de 1 1 inconvenient que la 

10 presence des sources ponctuelles augmente le dosage du 
rayonnement auquel le corps est soumis. De plus, pendant la 
collecte des donnees sur un cercle complet, les sources 
ponctuelles sont parfois situees a des positions ou le corps 
les cache du reseau de detecteurs. La traversee du corps 

15 peut attenuer fortement le rayonnement gamma, degradant 
I'aptitude a situer les sources ponctuelles. 

Une technique automatique de correction de mouve- 
ment dans un systeme angiographique est decrite dans un 
article de Alain Venot et V. LeClerk intitule "Automated 

20 Correction of Patient Motion and Gray Values Prior to Sub- 
traction in Digitized Angiography" , (Correction automatique 
du mouvement d'un patient et valeurs de gris avant soustrac- 
tion dans 1 1 angiographic numerisee), paru dans la revue IEEE 
Transactions of Medical Imaging, Volume Ml- 3, N° 4, Decembre 

25 1984, pages 179-186. Cette technique maximalise un critere 
de changement de signe deterministe par rapport a deux 
decalages en translation et une valeur constante. Lorsque 
les decalages et la valeur- constante sont tels que 1' image 
resultante, exempte de bruit, est proche de zero, tout bruit 

30 dans 1' image produit un decalage du signal, soit positif, 
soit negatif par rapport a zero. Au signal zero (bruit 
seulement), la probabilite du passage du signal du signe 
plus au signe rooins ou. vice versa est 0,5. Cela produit le 
nombre maximum de changements de signe.. Les valeurs non-op- 

35 timum du critere placent 1 1 image resultante plus loin du 
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zero et le bruit superpose donne une probability reduite de 
la production d'un changement de signe. Ainsi, la maximali- 
sation du changement de signe compense le mieux le mouvement 
du patient. 

Un autre article de Manbir Singh et al, intitule 
M A Digital Technique for Accurate Change Detection in 
Nuclear Medical Images - With Application to Myocardial 
Perfusion Studies Using Thallium-201 " (Technique numerique 
pour la detection precise de changements dans des images 
medicales nucleaires - avec applications aux etudes de 
perfusion du myocarde utilisant le Thallium-201 ) , paru dans 
la revue IEEE Transactions on Nuclear Science, Volume NS-26, 
N° 251, Pevrier 1979, pages 565-575, tente le cadrage 
d* images separees prises a des intervalles de, par exemple, 
une semaine, ou l*une des images est prise lors de la ten- 
sion creee par un exercice, et 1' autre en 1* absence d'une 
telle tension. Le procede decrit dans cet article necessite 
une interaction de la part de 1 ■utilisateur, contrairement 
aux techniques automatiques de la presente invention. 

En consequence, la presente invention a pour objet 
un moyen permettant de detecter le mouvement du corps entre 
des vues successives dans un systeme de balayage a vues 
multiples. 

La presente invention a pour autre objet un moyen 
et un procede pour corriger des donnees d 1 image afin d'eli- 
miner les erreurs dues au mouvement du corps pendant la 
collect e des donnees. 

La presente invention a encore pour objet un 
systeme d'imagerie d'un type produis£*nt une multitude de 
matrices de donnees X-Y representant des vues en projection 
d'un objet pour la reconstruction tomographique ulterieure 
de tranches axiales de 1' objet. Le mouvement detecte peut 
etre elimine en decalant ulterieurement de facon appropriee 
la donnee pour 1 'aligner avec une donnee anterieure, ou vice 
versa. Dans un mode de realisation de l a invention, toutes 
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les valeurs de brillance, paralleled aux axes X de deux vues 
se succedant, sont additionnees pour developper deux en- 
sembles a une dimension contenant une information sur le 
diagramme de brillance dans la direction Y. L'un des deux 

5 ensembles a une dimension est decale, et les deux ensembles 
a une dimension sont multiplies, element d ' image par element 
d' image, et les produits sont additionnes pour fournir une 
mesure de similarity entre les diagrammes de brillance dans 
les deux ensembles a une dimension. Des decalages repetes 

10 produisent une fonction de correlation liee a la position de 
decalage a laquelle ils se produisent. La position produi- 
sant un maximum dans la fonction de correlation est une 
bonne indication de la valeur du mouvement se produisant 
entre les deux vues. On peut utiliser le mouvement detecte 

15 pour signaler a l'operateur que la donnee est defectueuse, 
6u pour permettre un alignement automatique de toutes les 
donnees provenant d • un balayage. On peut executer un traite- 
ment identique en faisant la somme des donnees parallele a 
l'axe Y afin de detecter et de corriger le mouvement suivant 

20 I'axe X, dans certains cas limites. 

Pour detecter les decalages dans la direction Y 
(axial), on obtient un ensemble de valeurs a une dimension 
en additionnant, pour chaque valeur de Y, toutes les valeurs 
des elements d' image dans la direction X. Cet ensemble est 

25 ainsi une sequence de valeurs en fonction de Y. Dans le cas 
d' images • d 1 une projection ideale, les valeurs de cet en- 
semble ne sont pas fonction de 1' angle de vision par rapport 
a l'axe de rotation auquel la donnee est acquise. Ainsi, la 
donnee dans 1' ensemble Y est la meme, quel que soit 1' angle 

30 de vision. Dans le cas de la tomograpbie par emission, ou 
les effets de 1 'attenuation des photons sont importants, il 
y a des differences dans les valeurs de 1* ensemble Y en 
fonction de l 1 angle de vision. Cependant, comme la compa- 
raison s'effectue entre une image d'une vue et 1* image 

35 suivante, le procede de detection fonctionne egalement de 
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facon fiable dans ce cas. 

Pour detecter les decalages dans la direction X 
(tangentielle), un ensemble a une dimension de valeurs est 
obtenu en additionnant, pour chaque valeur de X, toutes les 
5 valeurs des elements d' image dans la direction Y. Cet en- 
semble est ainsi une suite de valeurs en fonction de X. Pour 
le cas d' images ideales de projection, les valeurs de cet 
ensemble changent generalement de facon considerable en 
fonction de I'angle de vision. Tout point unique dans l'es- 
10 pace-ob jet produit une variation sinusol'dale dans la direc- 
tion X en fonction de 1' angle de vision. Des images com- 
plexes produisent ainsi des ensembles de valeurs qui ne font 
pas l'objet d'une bonne correlation entre une vue et la vue 
suivante. Ainsi, en general, on ne peut utiliser le procede 
15 decrit pour detecter des decalages d'objet dans la direction 
x - Cependant, on a trouve que, s'agissant de la tomographie 
a emission, on rencontre souvent le cas ou il y a essentiel- 
lement une source distribute d'activite. Dans ce cas, le 
procede de correlation fonctionne bien en ce sens qu'il 
20 produit une valeur pour le decalage de donnees entre une vue 
et la vue suivante. On compare alors ce decalage au deplace- 
ment sinusoidal attendu et la difference fournit une mesure 
du patient dans la direction X. 

En bref, La presente invention fournit un systeme 
25 d'imagerie d'un type produisant une multitude de matrices 
X-Y representant des vues en projection d'un objet pour une 
reconstruction ulterieure utilisant des parties de la multi- 
tude de vues. Le systeme d'imagerie fait appel aux tech- 
niques de correlation pour determiner la valeur du mouvement 
30 suivant les axes X et Y. Le mouvement detecte peut etre 
elimine par un decalage ulterieur approprie de la donnee 
afin de I'aligner avec une donnee anterieure, ou vice versa. 
Dans un mode de realisation de 1* invention, toutes les 
valeurs de brillance paralleles aux axes X de deux vues sont 
35 additionnees afin de developper deux ensembles a une dimen- 
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sion contenant une information sur les diagrammes de bril- 
lance dans la direction Y. L'un des deux ensembles a une 
dimension est decale, et les deux ensembles a une dimension 
sont multiplies, element d' image par element d 1 image, et les 

5 produits sont additionnes pour fournir une mesure de simi- 
larity entre les diagrammes de brillance dans les ensembles 
a une dimension. Des decalages repetes produisent une fonc- 
tion de correlation, dorit le maximum est une bonne approxi- 
mation de la valeur du mouvement se produisant entre les 

10 deux vues. On peut utiliser le mouvement detecte pour si- 
gnaler a l'operateur que la donnee est defectueuse ou pour 
permettre un alignement automat ique de toutes les donnees 
provenant d'un balayage. On peut executer un processus 
identique suivant I'axe Y pour detecter et corriger tout 

15 mouvement suivant l'axe X, dans certains cas limites. 

Selon un mode de realisation de la presente inven- 
tion, on prevoit un dispositif dans un systeme d'imagerie 
pour corriger uh mouvement inter-vue, comprenantj un moyen 
pour recueillir des premiere et seconde vues d 1 un objet, ces 

20 premiere et seconde vues etant prises a partir de premier et 
second angles differents, chacune des premiere et seconde 
vues comportant au moins l'une des M rangees et N colonnes 
d'une donnee sur la brillance de 1' image, un moyen pour 
comparer un premier diagramme des donnees de brillance 

25 d 1 image en provenance de la premiere vue a un second dia- 
gramme correspondant des donnees de brillance d 1 image pro- 
venant de la seconde vue af in de produire un coefficient 
contenant une information sur une similarity entre eux, le 
moyen de comparaison comprenant un moyen pour appliquer un 

30 decalage a l'un des premier et second diagrammes pour une 
multitude de positions afin de produire une multitude des 
coefficients, un a chacune de la multitude de positions, et 
un moyen lie a la multitude de coefficient pour determiner 
une valeur du decalage donnant une similarite maximum. 

35 Selon une caracteristique de la presente inven- 
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tion, on prevoit un procede pour corriger le mouvement 
inter-vue dans un systeme d'imagerie, comprenant les etapes 
consistant a : recueillir des premiere et seconde vues d'un 
objet a partir de premier et second angles differents, 
5 respectivement, chacune des premiere et seconde vues compor- 
tant au moins 1'une de M rangees et N colonnes de donnees 
sur la brillance de 1' image, comparer un premier diagrarame 
de donnees de brillance de 1 1 image provenant de la premiere 
vue a un second diagramme correspondant des donnees sur la 
10 brillance de 1' image provenant de la seconde vue afin de 
Produire un coefficient contenant une information sur la 
similar ite entre eux, l'etape de comparaison comportant 
1 'application d'un decalage a l'un des premier et second 
diagrammes jusqu'a une multitude de positions afin de pro- 
15 duire une multitude de coefficients, un coefficient pour 
chacune des positions de la multitude, et determiner, sur la 
base de la multitude de coefficients, une valeur de decalage 
donnant la similarity maximum. 

La description qui va suivre se refere aux figures 
20 annexees qui representent respect ivement : 

Figure 1, un schema simplifie d'un systeme 
d'imagerie tomographique par emission ; 

figure 2, un schema auquel on se reportera lors de 
la description de la maniere avec laquelle on construit un 
25 sinogramme a partir de portions d'une multitude de vues dans 
le systeme d'imagerie tomographique de la figure 1 ; 

figures 3A-3C, des courbes simplifiees de bril- 
lance auxquelles on se reportera dans 1 ' explication de la 
technique de correlation employee dans la presente invention; 
30 figure 4, la fonction de correlation des courbes 

simplifiees de brillance des figures 3A-3C ; 

figure 5, un schema representant la technique de 
correlation utilisee pour obtenir la donnee de decalage dans 
deux vues se succedant ; 
35 figure 6, un ensemble de courbes representant 
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comment on peut analyser la donnee de correlation obtenue 
dans le dispositif de la figure 5 pour detecter le decalage 
devant etre applique a l'une des images ; 

figure 7, un schema fonctionnel sous forme de 

5 blocs d'un systeme d'imagerie employant un mode de realisa- 
tion de la presente invention. 

En figure 1, un patient 10 est supporte d'une 
maniere non representee au centre d'un trajet de balayage 12 
parcouru par un detecteur plan 14 de rayons gamma. Dans la 

10 preparation du balayage, le patient 10 regoit une injection 
d'une substance radioactive ayant une af finite etroite pour 
un organe 16 du corps, pour lequel on desire obtenir des 
donnees d' image. La substance radioactive se propage de 
maniere pref erentielle jusqu'a 1' organe 16, ce qui fait de 

15 cet organe une source radioactive, comme cela est indique 
par des rayons 18. Bien que les rayons 18 rayonnent dans 
toutes les directions, le detecteur 14, comme cela est 
classique, ne repond qu'aux rayons arrivant pratiquement 
perpendiculairement a son plan. 

20 Dans l'une de ses positions, representee en trait 

plein, le detecteur 14 detecte un ensemble rectangulaire M x 
N de.densites de rayonnement, chaque densite representant un 
element d'une image plane qu'on peut ulterieurement afficher 
sur un dispositif de visualisation approprie- M et N sont 

25 convent ionnellement egaux l'un a 1' autre et sont couramment 
de 64 x 64 ou .de 128 x 128 elements d' image. Chaque image 
%est appelee une vue. Les valeurs recueillies des elements 
d' image de la vue sont transmises a un" ordinateur 20 pour 
stockage et traitement ulterieurs, comme on 1 ' expliquera. 

30 Apres avoir complete une vue, le detecteur 14 

avance pas a pas sur le trajet de balayage 12 par increments 
d' angle predetermines jusqu'a une position illustree en 
tirets, ou il repond seulement au rayonnement tombant per- 
pendiculairement a son plan, represente par un rayon 18 en 

35 tirets. La vue obtenue dans cette position est de meme 
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transmise a 1'ordinateur 20. Ce processus s'effectue en 
deplacant le detecteur 14 suivant des angles egaux, sur un 
arc de 180 degres ou de 300 degres autour du patient 10, 
avec une prise de vue a chaque position, et en stockant des 
5 donnees dans 1'ordinateur 20* Les pas angulaires peuvent 
e tre, d' environ 5,6 degres de maniere a produire 64 vues 
pour un angle de 360 degres. II est classique de prendre 
autant de vues sur 360 degres qu'il y a d' elements damage 
dans une dimension d'une vue. 
10 Pour la description faite ulterieurement, on 

definit un axe Y comme etant un axe du detecteur 14 paral- 
lele a l'axe de rotation 22 de ce detecteur. On definit un 
axe X comme l'axe du detecteur 14 perpend iculaire a l'axe Y. 

En liaison maintenant avec la figure 2, on a 
15 represent! une vue 24 par un ensemble rectangulaire de 
rangees horizontales 26 et de colonnes verticales 28. Bien 
que les vues reelles 24 emploient classiquement 64 x 64 ou 
128 x 128 elements d' image, la reduction du nombre d 1 ele- 
ments d' image a un ensemble 8x8 diminue la complexity du 
20 dessin sans gener la comprehension du processus deer it. La 
troisieme rangee 26 des vues consecutives 24 forme le reste des 
rangees de 1* image du sinogramme 30. La donnee dans le 
sinogramme est alors utilisee pour reconstruire une image 
des tranches transaxiales en employant des techniques de 
25 reconstruction classiques telles que, par exemple, une 
projection par diffusion filtree. On peut construire de 
cette fagon autant d 1 images des tranches transaxiales qu'il 
y a de rangees 26. 

Comme la collecte des donnees dans les vues sue- 
30 cessives 24 est separee d'un laps de temps substantiel, le 
mouvement du corps dont on fait une image peut deplacer la 
donnee d'une vue 24 par rapport a celle des vues 24 voi- 
sines. II en resulte qu'une image transaxiale creee a partir 
de la donnee d'un sinogramme peut etre floue ou peut com- 
35 porter des elements ouvres d 1 image dus au mouvement qui 
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conduisent a un diagnostique manque ou faux. 

On a decouvert qu'on peut utiliser la donnee 
d' image dans des vues successives pour determiner la valeur 
du mouvement se produisant entre elles, et on peut employer 

5 1' information ainsi obtenue pour corriger la donnee et 
eliminer les effets d'un tel mouvement. 

Comme le detecteur plan 14 de rayons gamma (figure 
1) avance pas par pas entre une position et la position 
suivante, les vues adjacentes 24, bien que differentes, 

10 restent assez semblables quant au diagramme de brillance a 
celle qui precede. La similarite est suffisamment grande 
pour qu'il soit possible d'utiliser les techniques de corre- 
lation croisees entre vues adjacentes 24 pour detecter la 
presence d'un mouvement entre les deux vues 24, et deter- 

15 miner l'aroplitude de ce mouvement. 

En liaison avec la figure 3A, on a represente un 
diagramme simplifie de brillance pour une succession d 1 ele- 
ments d* image dans une colonne de vue N, la brillance etant 
limitee a. un maximum de quatre niveaux. Les second et 

20 septieme elements d 1 image ont une brillance egale a deux; le 
quatrieme element d' image une brillance egale a quatre, et 
tous les autres elements d' image une brillance egale a un. 
Une image reelle comporte beaucoup plus d 1 elements d' image 
et beaucoup plus de niveaux de brillance. 

25 En liaison maintenant egalement avec la figure 3B, 

on a illustre un diagramme <3e brillance correspondant dans 
la vue N+l f ou il ne s'est produit aucun decalage de 
1' image. On remarquera que le diagramme de brillance et 
d'obscurite de la figure 3B est semblable, mais non iden- 

30 tique, a celui de la figure 3A. Plus precisement, les 
pointes se produisent dans des elements d' image correspon- 
dant s. 

Le probleme, et une solution suggeree, sont illus- 
tres en figure 3C f laquelle represente la vue N+l dans 
35 laquelle la donnee a ete decalee de deux elements d' image 
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vers la droite. Dans le diagramme decale, la pointe dans 
l'element d 1 image 6, par exemple, est mal placee. Supposons 
que ces ensembles a une dimension representent les ensembles 
de valeurs qu'on obtient lorsqu'on fait I 1 addition dans la 
direction X des images des vues de projection. Alors, la 
dimension en abscisse, "element d" image" EI des figures 
3A-3C, correspond a la dimension Y decrite precedemment . 
Comme on l'a discute dans la description du sinogramme, la 
donnee gu'on utilise pour reconstruire une tranche trans- 
axial particuliere provient de la meme rangee Y dans toutes 
x es images des vues. Si la donnee de la figure 3A represent e 
la position Y correcte, alors celle de la figure 3C illustre 
la donnee provenant d'une image dans laquelle l'objet dont 
on fait 1* image s^st decale dans la direction Y par rapport 
a 1 1 emplacement correct. Si on essaie une reconstruction des 
tranches transaxiales en utilisant une telle donnee decalee 
par mouvement, la donnee choisie a partir de la figure 3C 
sera incorrecte pour cette tranche particuliere du corps. Au 
lieu de cela, on choisira la donnee provenant d'une autre 
tranche. 

La solution suggeree provient du fait que, bien 
que le diagramme de la figure 3C soit decale de deux ele- 
ments d 1 image vers la droite, sa forme n'est pas neanmoins 
tout a fait semblable a celle du diagramme de 1 1 image prece- 
dente en figure 3A. On a utilise cette similarity des formes 
pour decouvrir jusqu'ou on doit decaler le diagramme de la 
figure 3C pour trouver la meilleure adaptation avec le 
diagramme de la figure 3A. La valeur du decalage ainsi 
detect e peut alors etre appliquee a la donnee stockee nume- 
riquement pour aligner les images qui correspondent a la 
donnee decalee de la figure 3C avec celles correspondant a 
la donnee de la figure 3A, comme s'il ne s'etait produit 
aucun decalage. Cela elimine les elements de mouvement des 
images des tranches transaxiales reconstruites a partir de 
cette donnee. 
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Pour poursuivre la description, si -les positions 
de tous les elements d' image de la figure 3C sont decalees 
de plusieurs elements vers la droite ou vers la gauche et 
que la valeur de la brillance de chaque element d ' image dans 

5 le diagramme decale de la figure 3C est multipliee par la 
valeur de la brillance de chaque element d' image correspon- 
dant de la figure 3A, et qu'on additionne les produits de 
toutes les multiplications, la somme donne une mesure de la 
similar ite entre les deux diagrammes. Lorsque le decalage 

10 applique est tel que les diagrammes sont adaptes de maniere 
des plus etroites, on trouve un maximum pour la somme. Cela 
correspond au deplacement du diagramme decale de la figure 
3C jusgu'a la position occupee par le diagramme non decale 
de la figure 3B. 

15 En liaison maintenant egalement avec la figure 4, 

on a illustre les sommes des produits pour les diagrammes 
simplifies des figures 3A et 3C pour des decalages allant de 
-3 a +3 elements d' image. On remarquera que la somme trouvee 
avec un decalage de -2 elements d* image est sensiblement 

20 superieure a toutes les autres sommes. On devait s'y at- 
tendre car le diagramme de la figure 3C est decale de deux 
elements d • image dans le sens positif (vers la droite). On 
peut utiliser la pointe a un decalage de -2 tant pour indi- 
quer la presence d'un mouvement inter-vue que pour fournir 

25 une donnee pour la correction automatique de la donnee pour 
decaler le mouvement et restaurer la qualite de 1 ' image 
qu'on aurait obtenue s'il ne s'etait produit aucun mouvement. 

Les diagrammes. de brillance des figures 3A-3C sont 
beaucoup plus simples qu'une image reelle. De plus, la 

30 description simplifiee ignore la possibility de mouvement 
dans les deux dimensions. 

En liaison maintenant avec la figure 5, on a 
represents le procede pour 1' image plus complexe rencontree 
dans un monde reel. On conservera la simplification que 

35 permet la limitation des dimensions de la vue 24 a 8 x 8 
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elements d 1 image dans la description suivante. Les valeurs 
de brillance de tous les elements d 1 image dans chaque rangee 
26 de la vue N sont additionnees pour produire un ensemble 
32 de sommation a une dimension X, On execute une operation 

5 similaire sur chaque rangee 26 de la vue N+l. Les deux 
ensembles 32 de sommation X sont decales, puis multiplies, 
element d' image par element d' image, dans un correlateur Y, 
34, de la maniere decrite en liaison avec les figures 3A-3C, 
afin de detecter 1' apparition d'une pointe. Le decalage 

10 auquel on detecte une telle pointe indique la valeur du 
m Ouvement suivant l'axe Y se produisant entre les vues N et 
N+l. Un signal de decalage Y contenant cette information 
sort sur une ligne 36. En dehors de son utilite pour indi- 
quer la presence du mouvement, le signal de la ligne 36 peut 

15 egalement etre employe pour l'alignement d' elements d 1 image 
correspondants suivant l'axe Y dans les vues N et N+l comme 
si aucun mouvement ne s'etait produit. 

D'une facon similaire, les valeurs de brillance 
dans chaque colonne 28 de la vue 24 des vue N et N+l sont 

20 additionnees pour produire deux ensembles 38 de sommation Y 
qui sont decales, multiplies et additionnes dans un correla- 
teur X, 40, pour produire un signal sur une ligne 42 indi- 
quant la valeur du mouvement suivant l'axe X entre les vues 
N et N+l. 

25 Le nombre maximum des elements d' image sur les- 

quels on tente un decalage peut varier dans des applications 
differentes. Plus le nombre des elements d' image pour les- 
quels on entreprend la recherche d'une pointe est petit, 
plus le traitement est rapide. Avec une dimension d' element 

30 d' image d' environ 1/2 centimetre, on a decouvert qu'un 
mouvement depassant environ 10 elements d 1 image (de 1'ordre 
de 5 cm) se produit tree rarement. Cependant, dans le cas ou 
il se produit, il est tellement important qu'il est proba- 
blement preferable de recommencer l'imagerie. Ainsi, pour 

35 rendre la description plus concrete, et sans avoir pour 



2602602 



- 16 - 

intension de la limiter, on. emploie un decalage allant de 
-10 a +10 elements d' image dans la technique de correlation 
de la presente invention. 

II est peu probable que le mouvement reel entre 

5 vues soit egal a un nombre entier d'elements d' image. Plus 
precisement, le decalage peut etre, par * exemple, de 2,7 
elements d'image suivant l'axe X et de 5,3 elements d' image 
suivant l'axe Y. On a decouvert que la detection du mouve- 
ment jusqu'a une fraction d' element d' image est valable pour 

10 l'etablissement de la valeur exacte du decalage a appliquer 
a une vue. 

En liaison maintenant avec la figure 6, on a 
represent! un histogramme de correlation dans lequel la 
valeur 44 de correlation la plus elevee, representee en 

15 hachures horizontales; est flanquee de chaque cote par des 
valeurs de correlation contiguSs 46, indiquees par des 
hachures en diagonale. On emploie une courbe d'ajustement 48 
pour la valeur de correlation la plus grande, 44, et la 
valeur contiguS de correlation 46, dont la pointe 50 indique 

20 la position du decalage donnant la valeur de correlation 
maximum. 

En liaison maintenant avec la figure 7, un systeme 
d'imagerie, illustre dans ses grandes lignes par la refe- 
rence 52, emploie les techniques venant d'etre decrites.. Le 

25 systeme d'imagerie 52 comporte un dispositif 54 de collecte 
des donnees de base tel que, par exemple, un dectecteur plan 
14 de rayons gamma (figure 1) anime d'un mouvement de ba- 
layage pas par pas autour d'un patient (non represent! pour 
collecter un ensemble de M x N points de donnees a chaque 

30 position de balayage afin de produire une image d'une vue. 
Un calculateur de decalage X-Y, 56, qui peut fonctionner de 
la maniere indiquee en figures 2-6, produit des signaux de 
decalage X et Y pour application a un module de decalage 58. 
Bien qu'on puisse employer d'autres techniques dans le mode 

35 de realisation prefere de 1* invention, la donnee sur la 
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(N+l) eme vue est decalee de la maniere decrite, de maniere 
a 1 'aligner avec la N eme vue. Ainsi, apres la seconde vue, 
la N elne vue peut elle-meme etre une vue decalee. Lors du 
traitement de chaque vue, elle est ainsi alignee avec toutes 
les vues precedentes entre la premiere et 1 'avant-derniere. 

L a donnee corrigee ainsi produite est transmise a 
un module de stockage ou de visualisation 60 pour les func- 
tions indiquees. 

II n'est pas necessaire que les fonctions effec- 
tuees en figure 7, ou les fonctions de correction decrites 
anterieurement soient calculees pendant la collecte des 
donnees. Au contraire, et cela constitue le mode de realisa- 
tion prefere, la donnee due a toutes les vues est collectee 
et stockee. Le decalage est effectue separement. 11 peut 
etre avantageux, mais cela ne presente aucun inter et pour la 
Presente invention, de stocker en permanence la donnee non 
decalee d'origine. La presente invention n'empeche pas un 
tel stockage permanent et l'obtention de la donnee d* image, 
corrigee en mouvement, en provenance de ce stockage. 

Lorsque les images sont recueillies sur un arc de 
360 degres, la somme de tous les decalages image a image 
doit etre egale a zero. Des erreurs peuvent se produire dans 
les calculs des valeurs de decalage, rendant la somme des 
deplacements differente de zero. Les deplacements calcules 
dans le calculateur 56 peuvent etre totalises pour toutes 
les images pour un arc de 360 degres et tout ecart par 
rapport a zero etre stocke comme valeur d'erreur. On peut 
alors employer la valeur d'erreur dans le module de stockage 
ou de visualisation pour obtenir un decalage de polarisation 
applique a la donnee provenant de toutes les vues apres la 
premiere vue. La premiere vue reste non decalee et la donnee 
provenant des vues restantes est decalee conformement aux 
resultats du calcul dans le calculateur 56 plus ou moins le 
decalage de polarisation. 

Les techniques precedentes emploient * la donnee 
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relativement bruyante dans une vue pour comparaison a la 
donnee de la vue precedente afin de determiner la valeur de 
la correction necessaire. On a decouvert qu r une amelioration 
de la correction de 1* image peut etre obtenue en agissant 

5 sur la donnee collectee dans un processus a quatre etapes 
qui permet d'utiliser 1' ensemble complet des donnees collec- 
tees dans un balayage complet pour fournir des ensembles 
d* images sensiblement exemptes de bruit sur lesquelles on 
peut executer les techniques de _ detection de mouvement, 

10 eventuellement couplees aux techniques de. correction de 
mouvement . 

Au lieu de corriger le mouvement du patient sur 
chaque vue comme on l'a decrit precedemment, la technique 
perfectionnee emploie 1* ensemble complet des donnees d'ori- 

15 gine pour produire, par projection par transparence, la 
reconstruction des tranches transaxiales et leur nouvelle 
projection, une vue de reference etant liee a chaque vue 
dont le mouvement doit etre corrige. On peut employer des 
techniques de ponderation lors d'un tel traitement pour 

20 simuler les effets de 1 1 attenuation. Comme chaque vue de 
reference comprend des contributions de 1' ensemble entier 
des donnees plutot que d'une seule vue, le bruit dans la 
valeur de reference est sensiblement reduit. 

L 9 operation de correction de mouvement procede a 

25 la comparaison de la vue presente avec sa vue appropriee de 
reference pour determiner la valeur du mouvement de correc- 
tion necessaire. Plus precisement, si on agit sur la vue N 
pour determiner un decalage necessaire, 1 'autre vue corres- 
pondant a la vue N est employee pour la comparaison. Comme 

30 la vue de reference N est sensiblement exempte de bruit, et 
qu'on peut tenir compte des effets de 1 'attenuation dans son 
traitement, l'effet du bruit ou de 1 ' attenuation sur la 
determination du decalage est sensiblement reduit. La donnee 
stockee finale peut etre affichee de la maniere decrite. 

35 La technique perf ectionnee precedente fonctionne 
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sur la mime base de donnees des vues recueillies que la 
technique decrite en premier lieu. On trouve la difference 
dans les etapes supplementaires de traitement qui precedent 
1* operation de correction. Comme precedemment, 1 ' ensemble 

5 des donnees brutes collectees peut etre conserve pour emploi 
ulterieur. De fait, on envisage qu'un seul ensemble de 
donnees soit employe pour l'une ou I'autre des techniques ou 
pour les deux techniques sans necessiter une nouvelle pre- 
sence du patient. 

10 Dans certaines situations d' imager ie, 1' informa- 

tion desiree peut etre accompagnee d'une importante informa- 
tion de fond qui est constante, ou varie legerement, spatia- 
lement. II peut etre utile d'estimer une telle composante de 
fond du signal et de soustraire 1 ' estimation de la donnee 

15 d' image avant d'exeouter les calculs et deplacements de la 
figure 7. 
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RE VEND I CAT I ONS 

1. Dispositif pour corriger un mouvement inter- 
vues dans un systeme d'imagerie (52), caracterise en ce 
qu'il comprend : 

5 un moyen (54) pour recueillir des premiere et 

second vues (24) d'un objet ; 

les premiere et seconde vues etant prises a partir 
de premier et second angles dif ferents ; 

chacune des premiere et seconde vues comportant au 
10 moins 1 1 une de M rangees (26) et de N colonnes (28) - de 
donnees sur la brillance des images r 

un moyen pour comparer un premier diagramme de la 
donnee sur la brillance de 1 ' image provenant de la premiere 
vue a. un second diagramme correspondant de la donnee sur la 
15 brillance de 1 ' image provenant de la seconde vue afin de 
produire un coefficient contenant une information sur une 
simularite entre eux ? 

le moyen de comparaison comprenant un moyen (56) 
pour appliquer un decalage a l'un des premier et second 
20 diagrammes jusqu'a une multitude de positions afin de pro- 
duire une multitude des coefficients, un pour chacune des 
positions de la multitude de positions ; et 

un moyen lie a la multitude de coefficients pour 
determiner une valeur du decalage conduisant a la similarite 
25 maximum. 

2. Dispositif selon la revendication 1, caracte- 
rise en ce que le moyen de comparaison comprend s 

un premier moyen pour additionner une brillance de 
chaque element d* image dans l'une de chaque rangee (26) ou 
30 de chaque colonne (28) de la premiere vue (N) pour produire 
un premier ensemble de sommation a une dimension (32 ; 38); 

un second moyen pour additionner une brillance de 
chaque element d* image dans l'une de chaque colonne ou de 
chaque rangee dans la seconde vue (N+l) pour produire un 
35 second ensemble de sommation a une dimension (38 ; 32) ; 
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un moyen (34 ; 40) pour multiplier les valeurs de 
brillance d* elements correspondants dans les premier et 
second ensembles de sommation a une dimension afin d'obtenir 
une multitude de produits; et 
5 un moyen (34 ; 40) pour additionner la multitude 

de produits et obtenir ledit coefficient. 

3. Dispositif selon la revendication 1, caracte- 
rise en ce qu*il comprend en outre un moyen (58) repondant a 
la valeur du decalage pour decaler l'une des premiere et 

10 seconde vues suivant une quantite liee a la valeur de ce 
decalage. 

4. Dispositif selon la revendication 1, caracte- 
rise en ce que le moyen pour determiner une valeur du deca- 
lage comporte un moyen pour ajuster une courbe a une multi- 

15 tude desdits coefficients, une pointe de la courbe represen- 
tant ladite valeur. 

5. Dispositif selon la revendication 4, caracte- 
rise en ce que la multitude des coefficients comporte un 
coefficient maximum et des premier et second coefficients 

20 adjacents a ce maximum. 

6. Procede de correction d'un mouvement inter-vues 
dans un systeme d'imagerie, caracterise en ce qu'il comprend 
les etapes consistant a s 

recueillir des premiere et seconde vues d'un objet 
25 a partir de premier et second angles differents, respective- 
ment, chacune des premiere et seconde vues comportant au 
moins I'une de M rangees et N colonnes de donnees sur la 
brillance de 1 1 image ; 

comparer un premier diagramme des donnees sur la 
30 brillance de 1 1 image provenant de la premiere vue a un 
second diagramme correspondant des donnees sur la brillance 
de 1' image provenant de la seconde vue pour produire un 
coefficient contenant une information sur une similarite 
entre eux ; 

35 l'etape de comparaison comprenant 1 1 application 
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d'un decalage a l'un des premier et second diagrammes jus- 
qu'a une multitude de positions afin de produire une multi- 
tude des coefficients, un pour chacune des positions de la 
multitude de positions ; et 
5 determiner, sur la base de la multitude de coef- 

ficients, une valeur du decalage donnant la similarite 
maximum. 

7. Procede selon la revendicatioii 6, caracterise 
en ce que 1 1 etape de comparaison comprend les etapes consis- 
10 tant a : 

additionner une brillance de chaque element 
d' image de I'une de chaque rangee ou de chaque colonne dans 
la premiere vue pour produire un premier ensemble de sorama- 
tipn a une dimension ; 
15 additionner une brillance de chaque element 

1 image de 1 1 une de chaque rangee ou chaque colonne dans la 
second e vue pour produire un second ensemble de sommation a 
une dimension ; 

multiplier les valeurs de brillance d 1 elements 
20 correspondants dans les premier et second ensembles de 
sommation a une dimension afin d'obtenir une multitude de 
produits ; et 

additionner la multitude de produits afin d'obte- 
nir ledit coefficient. 
25 8. Procede selon la revendicatipn 6, caracterise 

en ce qu'il comprend en outre 1' etape consistant a decaler 
la premiere ou la seconde vue d'une quant ite liee a la 
valeur du decalage. 

9. Procede selon la reyendication 6, caracterise 
30 en ce que 1' etape de determination d'une valeur du decalage 

comporte 1' etape d'ajustement d'une courbe a une multitude 
desdits coefficients, une pointe de la courbe representant 
ladite valeur* 

10. Procede selon la revendication 9, caracterise 
35 en ce que la multitude desdits coefficients comporte un 
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maximum et des premier et second coefficients adjacents a ce 
maximum. 

11. Dispositif pour la correction d'un mouvement 
inter-vues dans un systeme d'imagerie, caracterise en ce 
qu * i 1 comprend r 

u n moyen pour recueillir et stocker une multitude 
de vues originales d'un objet prises dans un cercle I'en- 
tourant ; 

chacune des vues originales comportant M rangees 
et N colonnes de donnees sur la brillance de 1* image ; 

un moyen pour proceder a une projection par trans- 
parence de la multitude de vues originales pour produire un 
ensemble de donnees a trois dimensions representant 1' objet ; 

un moyen pour proceder a une projection par trans- 
parence de 1' ensemble de donnees a trois dimensions afin de 
produire une multitude de vues de reference ; 

chacune des vues de reference comportant M rangees 
et N colonnes correspondant aux vues originales, 

chacune des vues originales ayant une vue corres- 
pondante de reference ; 

un moyen pour comparer un premier diagramme de la 
brillance de 1 ' image provenant d'une premiere vue originale 
a un second diagramme de la brillance de 1 1 image provenant 
de sa vue de reference pour produire un coefficient conte- 
nant une information sur une similarite entre eux ; 

le moyen de comparaison comportant un moyen pour 
appliquer un decalage au premier diagramme jusqu'a une 
multitude de positions afin de produire une multitude des- 
dits coefficients, un coefficient pour chaque position de la 
multitude desdites positions ? et 

un moyen lie a la multitude de coefficients pour 
determiner une valeur du decalage donnant la similarite 
maximum. 

12. Precede pour corriger le mouvement inter-vues 
dans un systeme d'imagerie, caracterise en ce qu'il comprend 
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les etapes consistant a : 

recueillir et stocker une multitude de vues ori- 
ginales d'un objet prises suivant un cercle 1 ' entouran-t ; 

chacune des vues originales comportant M rangees 
5 et N colonnes de donnees sur la brillance de 1' image ; 

proceder a une projection par transparence de la 
multitude de vues originales pour produire un ensemble de 
donnees a trois dimensions " represantant 1' objet ; 

proceder a une projection par transparence de 
10 l'ensmble de donnees a trois dimensions pour produire une 
multitude de vues de reference ; 

chacune des vues de reference comportant M rangees 
et N colonnes correspondant aux vues originales ; 

chacune des vues originales ayant une vue corres- 
15 pondante de reference ; 

comparer un premier diagramme de la brillance de 
1' image provenant d 1 une vue originale a un second diagramme 
de brillance de 1 1 image provenant de sa vue de reference 
pour produire un coefficient contenant une information sur 
20 une similarity entre eux ; 

appliquer un decalage au premier diagramme jusqu'a 
une multitude de positions afin de produire une multitude 
desdits coefficients, un coefficient pour chaque position de 
la multitude de positions ; et 
25 determiner a partir de la multitude de coeffi- 

cients une valeur du decalage indiquant la similarite .maxi- 
mum . 
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